Xkmeans = [gauss_sample([0;0], [1,0; 0,1], 100), gauss_sample([5;0], [1,0;0,1], 100), gauss_sample([0;5], [1,0;0,1],100)]';
clusters = k_means(Xkmeans, 3);

figure
scatter(Xkmeans(:,1), Xkmeans(:,2), 'x', 'k');

figure
hold on
scatter(Xkmeans(clusters == 1,1), Xkmeans(clusters == 1,2), 'x', 'r');
scatter(Xkmeans(clusters == 2,1), Xkmeans(clusters == 2,2), 'x', 'g');
scatter(Xkmeans(clusters == 3,1), Xkmeans(clusters == 3,2), 'x', 'b');